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ABSTRACT 

We report the discovery of a new 21-cm Hi absorption system using commissioning 
data from the Boolardy Engineering Test Array of the Australian Square Kilometre Ar¬ 
ray Pathfinder (ASKAP). Using the 711.5-1015.5 MHz band of ASKAP we were able 
to conduct a blind search for the 21-cm line in a continuous redshift range between 
z = 0.4 and 1.0, which has, until now, remained largely unexplored. The absorption line 
is detected at z = 0.44 towards the GHz-peaked spectrum radio source PKS B1740—517 
and demonstrates ASKAP’s excellent capability for performing a future wide-held sur¬ 
vey for Hi absorption at these redshifts. Optical spectroscopy and imaging using the 
Gemini-South telescope indicates that the HI gas is intrinsic to the host galaxy of the 
radio source. The narrow [O III] emission lines show clear double-peaked structure, in¬ 
dicating either large-scale outflow or rotation of the ionized gas. Archival data from the 
XMM-Newton satellite exhibit an absorbed X-ray spectrum that is consistent with a high 
column density obscuring medium around the active galactic nucleus. The HI absorption 
profile is complex, with four distinct components ranging in width from 5 to 300 km s -1 
and fractional depths from 0.2 to 20 per cent. In addition to systemic HI gas, in a circum- 
nuclear disc or ring structure aligned with the radio jet, we find evidence for a possible 
broad outflow of neutral gas moving at a radial velocity of v ~ 300 km s _1 . We infer that 
the expanding young radio source (t age « 2500 yr) is cocooned within a dense medium 
and may be driving circumnuclear neutral gas in an outflow of ^ 1 M© yr -1 . 

Key words: methods: data analysis - ISM: jets and outflows - galaxies: active - galax¬ 
ies: ISM - radio lines: galaxies. 
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1 INTRODUCTION 


Over the past 60 yr the 21-cm hyperfine transition of atomic hy¬ 
drogen (H I) has been used by astronomers to measure the distri¬ 
bution and kinematics of the neutral interstellar medium (ISM) 
in galaxies and, importantly, trace the available fuel for fu¬ 
ture star-formation. Wide-held surveys have systematically mea- 


20011 Meyer et al. 2004, Walter et al. 2008 

Koribalski 

2010 

Oosterloo et al. 2010, Haynes et ak|2011, 

eald et al. 

2011 

Serra et al. 2012, Wang et al. 2013), while at greater distances 


individual Hi-rich galaxies have been detected out to z ~ 0.3 
(e.g. Catinella et al.|2008]|Verheijen et al.|2010||Freudling et al 


2011). Statistical techniques, such as spectral stacking (e.g. Lah 


et al.|2 009, Del haize et al.|2013"l|Rhee et al. 2013 Gere b et al. 


2014| ) and intensity mapping (e.g. Chang et al. 2010 , Masui et al. 
2013) have also provided constraints on the cosmological Hi 


mass density at z < 0.4 and z < 1.0 respectively. 

At higher redshifts it becomes increasingly difficult for ex¬ 
isting radio telescopes to detect the faint 21-cm emission from 
individual galaxies and we must look to absorption against back¬ 
ground sources, either at 21-c m (e.g.|Carilli et al.|1998]|Curran| 


|et al. |2006[[2008)|20 11 [ |20 1 3b[ IKanekar et al.|2009t oT the Ly a 

line at 1216 A (e. g.|Peroux et al.|2003 Rao et al.|2006 Noter- 


|daeme et al.|2012||Zafar et al.|2013| ), to understand the evolving 
HI content of the distant Universe. Deficiency of spatial infor¬ 
mation is the significant drawback of using absorption to map 
the neutral ISM, yet it still provides one of the few practical 
methods of directly detecting cold HI gas in individual galaxies 
at an epoch of the Universe (0.5 < z < 2) that has, until now, 
remained largely unexplored (e.g. [Lagos et al.|20T4| . 

Aside from intervening galaxies, HI absorption is also a 
direct tracer of the interaction between radio source and host 
galaxy. Recent 21-cm surveys of radio-loud active galactic nu¬ 
clei (AGN) have revealed the varied kinematic signatures of ro¬ 
tating discs, infalling and outflowing gas, and individual offset 


clouds (e.g. van Gorkom et al. 1989 

Morganti et al. 2001 Ver- 

meulen et al. 2003, Morganti et al. 

2005, |Gupta et al. 

2006 

Emonts et al. 2010, Chandola et al. 2011, Allison et al. 

2012, 

2014, Gereb et al. 2015). Amongst these heterogeneous samples 


we often find that the most compact sources have the highest 
rate of detection, the result of either intrinsically higher column 
densities of neutral gas ( jPihl strom et al.|2003), the relative ori¬ 
entation between source and absorber ( [Curran et al.|2013a| ) or, 
quite possibly, a combination of both (Orienti et al. 2006, Gereb 


et al. 2015). The population of compact steep-spectrum (CSS; 
d < 15kpc) and gigahertz peaked-spectrum (GPS; d < lkpc) 
radio sources are particularly interesting targets since they are 
thought t o constitute recently triggered radio AGN (t 
- 10 6 yr; 


Fanti et al. 1995 


Readhead et al. 1996 Owsianik 


age_ 


10 4 


|& Conway 1 1998]) . Higher detection yields of HI absorption in 
these sources may be indicative of the dense and dusty environ¬ 
ments in which they are born. 

By studying the kinematics of the HI gas in young radio 
galaxies, along with the molecular and ionized gas, we can un¬ 
derstand the processes by which these radio jets are triggered 
and their subsequent feedback on the neutral ISM. |Vermeulen] 
|et al.| ( |2003|) and |Morganti et al.| ( |2005| ) found examples of 
blue-shifted HI absorption components at large negative veloc¬ 
ities (~ —1000 km s -1 ) with respect to galaxies hosting com¬ 
pact radio sources, symptomatic of jet-driven outflows of neu¬ 
tral gas. More recently, in-depth case studies of individual ra¬ 
dio AGN, using high spatial resolution 21-cm and multi wave¬ 
length data, have provided direct evidence for the existence of 
~ 10 M 0 yr -1 jet-driven outflows of Hi, co-located with out¬ 
flows of ionized and molecular gas e.g. 1504+377 ( jKanekar ~&\ 


Chengalur 

2008 

hony et al. 

2013 


and IC 5063 ( Tadhunter et al.|20l4) . 


The advent of precursor and pathfinder telescopes to the 
Square Kilometre Array (SKA), in particular the Australian 
SKA Pathfinder (ASKAP; | Johnston et al.||20071 |Deboer et al.| 

|2009[|Schinckel et al.|2012| ), the South African MeerKAT tele¬ 

scope ( |Booth et al.||2009| and the We sterbork Aperture Tile 
in Focus (APERTIF; |Oosterloo et al.||2009| >, will enable as¬ 
tronomers to carry out radio-selected surveys for HI absorption 
over most of the sky. Crucially, such surveys are made possi¬ 
ble by significant improvements in bandwidth, field-of-view, and 
new observatories with reduced terrestrial radio frequency inter¬ 
ference (RFI). 

Here we report the first result from a search for redshifted 
absorption using the ASKAP Boolardy Engineering Test Ar¬ 
ray (BETA; [Hotan et al.||2014| ), a six-antenna prototype de¬ 
signed to demonstrate the feasibility and science capability of 
phased-array feeds (PAFs) based on a planar connected “che- 
querboard” array ( [Hay & O’Sullivan|20Q8| ). The backend digital 
beamformer weights the signal from 188 independent PAF re¬ 
ceptors at the focal plane of each antenna to electronically form 
up to nine simultaneous primary beams within an area of ap¬ 
proximately 30 deg 2 . This impressive flexibility in field-of-view 
is matched by the telescope’s spectral capability; the fine filter- 
bank generates 16416 channels over 304 MHz of bandwidth, at 
observing frequencies between 0.7 and 1.8 GHz, achieving an 
average 21-cm line resolution of 5.5kms _1 and maximum Hi 
redshift of z = 1.0. We are currently using BETA to search 
for HI absorption towards the brightest and most compact radio 
sources in the southern sky. We will continue to modify our tar¬ 
get sample to fainter sources as more antennas are added to the 
ASKAP array, culminating in the First Large Absorption Survey 
in Hi (FLASH). 

Throughout this paper we adopt a flat A cold dark matter 
(ACDM) cosmology with Ho = 70kms -1 Mpc -1 , Qm = 0.3 
and £2 a = 0.7. Radial velocities and redshifts have been corrected 
for the solar barycentric standard of rest frame. Uncertainties are 
given as 1 a intervals unless otherwise stated. 


2 HI ABSORPTION WITH BETA 

2.1 A pilot survey of bright, compact radio sources 

As part of the early commissioning and science demonstration 
phase of ASKAP, we are using BETA to carry out a pilot survey 
of Hi absorption at 0.4 < z < 1.0 towards the brightest and 
most compact radio sources in the southern sky. The sensitivity 
of the telescope is substantially lower than will be achieved with 
the full ASKAfQand so we select targets that optimize the op¬ 
tical depth sensitivity, through both a high signal-to-noise ratio 
(S/N) continuum background and high expected fraction of ra¬ 
dio flux obscured by foreground neutral gas. The following is a 
list of criteria from which we select our targets. 

(i) Visible at declinations of S < +20°. 

(ii) A total flux density greater than 1 Jy in existing all-sky 
catalogues at similar frequencies to this band, which include the 
National Radio Astronomy Observatory Very Large Array Sky 
Survey (NVSS; [Condon et al.|1998|> , the Sydney University Mo- 
longlo Sky Survey (SUMSS; |Mauch et al.|2003| and the second 
epoch Molonglo Galactic Plane Survey (MGPS-2; Murphy et al. 
[2507] ). 


1 The full array is expected to have 36 antennas, each fitted with a PAF 
that can simultaneously form up to 36 beams. 
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Figure 1. The SED at radio wavelengths of GPS radio source 
PKS B1740—517. The frequency axis is given in the observer rest- 
frame. The dashed line denotes a best-fitting broken power-law model 
to the data (Moffet|1975) . Using standard least squares minimization we 
obtain a break frequency of 1 GHz and optically thick and thin spec¬ 
tral indices of a = 0.83 and —0.74, respectively. References for the 


data: Jacll -| Jacob s et al.H2011 1 ; Lar81 - |Large et al.||1981|; Cam94 - 
|Campbell - Wilson & Hunsteadfi 1994) ; Mau03 - |Mauch e tal. 1(2003) ; 


Tin03 - Tingay et al. <2003); Wal75 - (Wall et al. 
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; McC12 -[Me: 

Connell et al. < 2012); Gregory et al. (1994) - Gre94; Hea07 - Healey 

et al. l 

20071; Ric06 - 

Ricci et al. ( 2006); Lop07 -|Lopez-Caniego et al. 

<2007 

Plall -Plane 

k Collaboration VII (20111. 


^ i i i i i i | 

A 


^iiiii 


★ 

Jacll 

* 

Wri94 

• 

Lar81 

◄ 

McC12 

T 

Cam94 

► 

Gre94 

A 

Mau03 

X 

Hea07 

• 

Tin03 

# 

Ric06 

□ 

Wal75 

+ 

Lop07 

♦ 

Mas08 

j _i_i_ 

# 

_i _ i i i i 

Plall 

J_L_ 


t 


*\ 


_i_i i i i I_i_i_ 


(iii) Either z ^ 0.4 or no known redshift according to the 
NASA/IPAC Extragalactic Database (NED^ thus excluding 
those sources that are known to be located in front of the vol¬ 
ume probed by 0.4 < z < 1.0. 

(iv) A significant fraction of the total flux density is dis¬ 
tributed on the compact scales measured by very-long-baseline- 
interferometry (VLBI), based on those sources listed in the 
VLB I radio fundamental catalogu^] 

We present here our first new detection of 21-cm absorption us¬ 
ing BETA as a demonstration of the capability of the ASKAP 
telescope to perform radio-selected surveys for redshifted HI. 
In future work we will present the results of further observations 
from this pilot survey. 


2.2 PKSB1740—517 

PKS B1740—517 was one of the first sources selected for BETA 
observations based on the selection criteria outlined above. It 
is an archetypal GPS radio source ( [Randall et aH2011| ), peak¬ 
ing at an observed frequency of 1 GHz, with optically thick and 
thin spectral indice^j of a — 0.83 and —0.74, respectively 
(see Fig.[lJ. The relative compactness and high flux density 
($843 = 8.15=b0.24 Jy; |Mauch et al.|2003) of this source make it 
an ideal target for HI absorption. Flux monitoring by King et al. 
over ~700 d, using the 26-m antenna at the Mt. Pleasant Obser¬ 
vatory (King 1994; Jauncey et al. 2003), indicated that the source 
is weakly varying at 2.3 GHz (by approximately 20 per cent), 
with no variation seen at 8.4 GHz over this time range. They 


2 http://ned.ipac.caltech. edu/ 

3 http://astrogeo.org/rfc jPetrov] 2013] and references 
therein). 

4 We adopt the convention that the spectral index a is defined by S' oc 
, where S is the flux density and v is the frequency. 


demonstrated, by fitting to visibility data from three baselines 
in the network of Australian VLBI antennas (Parkes - Hobart, 
Parkes - Perth and Parkes - Alice Springs), that the 2.3 GHz con¬ 
tinuum emission is adequately modelled by two compact Gaus¬ 
sian components of angular extent less than 10 mas and sepa¬ 
rated by 52 mas (Fig.|2]. While there is some residual flux on the 
longest baselines, which might be indicative of more complex 
compact structure, the total flux density at 2.3 GHz is apparently 
accounted for by the modelled compact components. The spec¬ 
tral behaviour, VLBI-scale structure and low variability are all 
evidence that this radio source is intrinsically young and/or con¬ 
fined by its environment (e.g. |Q’Dea|1998) . 

At a Galactic latitude of b = —11? 5 and longitude of 
l = 340?2, PKSB1740—517 is seen through a densely popu¬ 
lated foreground (Fig.[2}, which is likely the reason for the lack 
of optical information for this source. Despite this, |di Serego") 
Alighieri et al. (1994) used imaging and spectroscopy with the 
ESO 3.6-m telescope to secure an optical identification of the 
host galaxy, noting that it has a particularly red continuum. Un¬ 
fortunately they could not detect any strong optical emission 
lines and hence secure a spectroscopic redshift. Approximate 
indicators of the redshift have been determined by photomet¬ 
ric means: Wall & Peacock (1985) give a rough estimate of 
z = 0.347 based on their empirically derived U-band magni¬ 
tude relationship (but noting a poor optical identification with 
the radio source), while [Burgess & Hunstead] |2006| estimate 
that z = 0.63 based on the 77-band magnitude of 20.8. At longer 
wavelengths, the infrared colours [3.4 - 4.6 fi m] = 1.06 and [4.6 
- 12 /i m] =2.85, from the Wide-field Infrared Survey Explorer 
(WISE; [Wright et aL||2010|>, are consistent with a quasi-stellar 
object at z ~ 0.5 (fig. 1 of |Blain et al.|2013) . These photometric 
indicators of the redshift suggested that PKS B1740—517 would 
provide a suitable background target against which we might de¬ 
tect absorption in the range 0.4 < z < 1.0 with BETA. 


2.3 Observations 

On three dates - 2014 June 24, August 03 and September 01 
- we used the 711.5-1015.5 MHz band of the ASKAP BETA 
prototype ( |Hotan et al.|2014| ) to search for 21-cm absorption to¬ 
wards PKSB1740—517 in the redshift range 0.4 < z < 1.0. 
In Table[l]we summarize the parameters for each of these ob¬ 
servations. Fig.[3] shows the nominal footprint of the nine PAF 
beams for our observation (a symmetric diamond pattern), cen¬ 
tred on RA(J2000) = 17 h 44 m 25?45 and Dec.(J2000) = 
—51°44 / 43 / /8. Additional observations of a calibrator source 
(PKS B1934—638) were carried out in short 5 - 15 min scans 
at the centre of each PAF beam. Using weights for the PAF ele¬ 
ments that maximize the S/N in the required direction (see Hotan 
[eTak] 20l4), beams were generated with nominal half-power 
widths that decrease from 2?07 to 1?45 across the whole band¬ 
pass. The six 12 m antennas of BETA are arranged in an array 
that is elongated by a factor of approximately 2 in the east—west 
direction, with baselines in the range 37 - 916 m. Therefore the 
data are sensitive to angular scales in the range 1-39 arcmin 
across the band. The fine filterbanks in the beamformer gener¬ 
ate 16416 independent spectral channels separated by approxi¬ 
mately 18.5kHz, equivalent to Hi velocities in the range 5.5 - 
7.8kms -1 . 


© 2015 The Authors, MNRAS 000,[TJ{T8] 
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Figur e 2. Left: an R -band optical image centred on PKS B1740—517 (300 s exposure with the 33 m An glo-A ustralian Telescope; |Burgess & Hunstead] 
|2006) , overlaid with SUMSS 843 MHz radio contours (5.5, 6.0, 6.5,7.0 and 7.5 Jy beam -1 ; Mauch et al .|2003) . The open red square denotes the position 
of the radio source - RA(J2000) = 17 h 44 m 25?450 and Dec.(J2000) = —51°44 / 43 / /79 - given by the second realization of the International 
Celestial Reference Frame (ICRF2; |Fey et al.|2009| >. The faint background streaking in the optical image is due to scattered light from /j, Arae, a fifth 
magnitude star at an angular separation of 5.9arcmin to the south-west. Right: a high spatial resolution 2.3 GHz source model of PKS B1740—517, 
obtained by |King|{l994) from VLBI observations using three antenna baselines, taken on 1989 December 06. The positions are given relative to the 
centroid of the brighter component. The contour levels are at 1, 5, 10, 20, 35, 50, 65, 80 and 95 percent of the 2.67 Jy beam -1 peak flux density. The 
model components have been convolved by the synthesized beam given by the shaded ellipse in the bottom left-hand corner. The source is well modelled 
by two hotspot components separated by 52 mas in a flux ratio of 3.9 to 1. 


2.4 Data analysis 

We used the CAS£^] package (McMullin e t al. 2007 ) to sepa¬ 
rate the data into each of the nine beams and flag the auto¬ 
correlations and large amplitude values due to digital glitches. 
The PAF beams are formed digitally by applying weights to el¬ 
emental receptors, which in the case of BETA are constant in 
a repeating pattern of 4 and 5 MHz intervals. There are some 
sharp bandpass effects at the edges of these intervals, but within 
them we find that the bandpass is smooth for most of the band at 
a level consistent with the noise in our data. This is further cor¬ 
roborated by the absence of any false positive detection; see Sec- 
tion |3.3| We therefore pursued a calibration strategy that splits 
the band, at full spectral resolution, into 64 chunks based on the 
exact pattern of 4 and 5 MHz intervals. Separately we generated 
a single data set of coarse 32 x 9.5 MHz channels across the 
band by averaging over 513 of the fine channels. The 18.5-kHz 
fine channels are used to search for Hi absorption, while the 

9.5 MHz coarse channels produce high S/N images that are used 
for self-calibration. 

We performed calibration, imaging and further gagging of 
the data using tasks from the MI RIAE[j package (Sault et al. 


|1995| ). Those 18.5-kHz channels that were corrupted by narrow- 
band RFI were flagged using a MIR I AD implementation of the 
SumThreshold method developed by |Offringa et al.| ( |2010| ), 
which typically resulted in data loss of a few per cent. The ma¬ 
jority of this low-level RFI is caused by single-channel spikes, 
commonly referred to as birdies, which are generated by the 
cooling system flow regulator valves associated with the first 
generation of phased-array feeds. Further corrupted data, result¬ 
ing from problems with individual antennas, antenna baselines 
or the backend correlator, were also manually removed at this 


5 http://casa.nrao.edu 

6 http://www.atnf.csiro.au/computing/software/ 
miriad/ 


stage. In particular, there was an intermittent problem with the 
time synchronisation between individual correlator cards that re¬ 
sulted in 50 per cent of the spectral band being flagged for the 
data obtained on 2014 June 24 and September 01. 

The flux density scale and rough gain corrections for 
each antenna were calculated using short scans (between 5 and 
15 min) of PKS B1934—638 in each PAF beam, based on the 
model of |Reynolds \\ 1994] ). These solutions were then transferred 
to the corresponding data in the PKS B1740—517 field. We 
achieved further correction for the time-dependent gains by per¬ 
forming multiple iterations of self-calibration on the 9.5 MHz 
coarse-channel data, which span the whole band and thereby 
provide optimal S/N. Initially the SUMSS 843 MHz radio cat¬ 
alogue ( [Mauch et al.|2003] ) was used to construct a sky model 
for each PAF beam, followed by models generated from the data 
using a multi-frequency synthesis variant of the CLEAN algo¬ 
rithm |Hogbom||l974| . At each iteration the MIRIAD SELF- 
CAL task was used to solve for the time-dependent gains, and 
we simply transferred these corrections to each of the corre¬ 
sponding 64 chunks of 18.5 kHz fine-channel data. 

Precise subtraction of the continuum flux density was 
achieved by generating continuum images for each of the 64 
chunks of fine-channel data and then subtracting the CLEAN 
component models (using UVMODEL) from the visibilities. 
Removal of any residual flux density, including curvature of the 
instrumental bandpass towards low frequencies, was performed 
by fitting a second-order polynomial to each chunk and sub¬ 
tracting using UYLIN. Channels containing any detected ab¬ 
sorption features were subsequently excluded in iterative refine¬ 
ment of the above continuum subtraction procedure. We im¬ 
age the continuum-subtracted data in the standard way, forming 
cubes that are equal to the nominal full width at half-maximum 
(FWHM) of the beams at the middle of the band. Spectra were 
then extracted at the positions of identified continuum target 
sources within each beam centre and converted to units of frac¬ 
tional absorption through dividing by the continuum. 


© 2015 The Authors, MNRAS 000,[l]{T8] 
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RA (J2000) 

Figure 3. A total intensity image of the PKSB1740—517 field observed with BETA (grey-scale; data obtained from a commissioning observation 
carried out on 2014 September 01). The small blue circles denote the positions of sources brighter than 100 mJy within l°of each beam centre, listed in 
the 843 MHz SUMSS jMauch et al.|2003) and MGPS-2 {Murphy et al.[2007| i catalogues. The large black circles denote the nominal beam width at half 
power at 863.5 MHz (1.028A/Ddi s h ~ 1.7°) for each of the nine beams. 


Table 1. Summary of our 21-cm observations using BETA, where column 1 gives the date of observation; column 2 and 3 the modified Julian 
date (MJD) start and end times; column 4 the number of PAF beams used and their footprint; column 5 the ASKAP antennas used; column 
6 the on-source integration time; columns 7-9 properties of the restoring beam, which include the FWHM of the major and minor axes, and 
their position angle, respectively; column 10 the median per-channel noise in the centre beam, across the 304 MHz band. Note that differences 
in the spectral noise between these observations are consistent with the integration times and number of antennas. 


Date 

MJD 

PAF beams® 

Antennas ' b 

fint 

$max 

^min 

<t> 

^chan 


Start 

End 



(h) 

(arcsec) 

n 

(mJy beam -1 ) 

2014 June 24 

56832.442 

56832.921 

1 (centre) 

1,8,9, 15 

11.5 

90 

70 

+70 

23 

2014 August 03 

56872.649 

56872.774 

9 (diamond) 

1,3,6, 8, 15 

3.0 

340 

100 

-40 

37 

2014 September 01 

56901.224 

56901.723 

9 (diamond) 

1,6, 8,9, 15 

11.9 

100 

80 

-90 

18 


a See Fig.[3]for details of the 9-beam diamond configuration. The centre configuration refers to a single beam at the pointing centre. 
b See fig. 2 of Hotan et al. (2014) for details of the BETA antenna positions. 


2.5 Combining spectra from multiple PAF beams 

Some sources in the field may be common to multiple PAF 
beams (see Fig(3j and so we construct a single spectrum by av¬ 
eraging over spectra from the individual beams. Since our ar¬ 
rangement of the PAF beams here is such that they do overlap, 
the noise in each spectrum is not independent and so the opti¬ 
mal S/N is given by weighting each spectrum by the inverse of 
the noise covariance matrix. In the case of the diamond footprint 
used here, the nominal noise correlation measured between adja¬ 
cent beams (separated by 1? 2 between their centres) ranges from 
approximately 30 to lOpercent across the 711.5 - 1015.5 MHz 


band. The optimal averaged spectrum is then given by 

& = &S [ W Cbeam &] (1) 

and 

O’! = [w t Cbe am w]“f ( 2 ) 

where Cbeam is the noise covariance matrix between beams, S 
is the set of fractional absorption data in each beam and w is an 
accompanying set of weights, which in this case are unity. The 
elements of Cbeam are estimated empirically by measuring the 
noise covariance between pairs of beams on a per-pixel basis. 

In Fig. [4] we show the resulting full 711.5-1015.5 MHz 
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Figure 4. The 711.5-1015.5 MHz BETA spectrum towards PKS B1740—517, averaged over all three observing epochs. For visual clarity the data have 
been binned from the native spectral resolution of 18.5 to 100 kHz. The barycentric corrected observed frequency is shown on the lower abscissa, and 
the upper-abscissa denotes the corresponding HI redshift. The data (black line) denote the change in flux density as a fraction of the continuum and the 
grey region gives the corresponding rms spectral noise multiplied by a factor of 5. An absorption line is visible in the spectrum at z'bary = 985.5 MHz, 
equal to a Hi redshift of z = 0.4413. In Fig.[5]we show this absorption line at full spectral resolution. An increase in the noise is seen in those parts 
of the band where 50 per cent of the correlator was unavailable for two of our three observations. Other noise features are the result of flagging and 
occasional failures of individual correlator cards. 
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Table 2. A summary of model parameters derived from fitting the Hi absorption line seen towards PKS B1740—517 for each epoch and 
the average spectrum. Column 1 gives the observation epoch; column 2 the Gaussian component corresponding to that shown in Fig.[5] 
column 3 the component redshift; column 4 the component rest-frame FWHM; column 5 the peak component depth as a fraction of the 
continuum flux density; column 6 the improvement in Bayes evidence for each component relative to the noise-only model (see |Allison| 
|et al.|2014] and references therein); column 7 the reduced chi-squared statistic for the best-fitting model parameters after the component 
is introduced into the model. Intervals of 1 a are given for the measured uncertainties, derived from simultaneously fitting all model 
Gaussian components and taking into account the gain response of the spectral channels (see text for details). 


Epoch 

ID 

-^bary 

AU 50 

(AS)/ »Scont)peak 

Aln(Z) 

X 2 /d.o.f. 




(kms -1 ) 

(per cent) 



2014 June 24 

1 

0 44129264 +0 ' 00000058 
u.^rzyzo^_ 0 00000060 

c 1 ir+0.20 
°- iO -0.21 

-20.201°;®® 

893.47 ±0.08 

3.29 


2 

fl 441 999Q+0-0000023 
U. z i z ilZZZo_ 0 0000022 

6-8+9 

—4 21 4-0 ' 61 
4.ZD _ 2 46 

99.59 ±0.14 

1.96 


3 

fl 441 £1 7+0-000016 
U.441ol7_ 0 000015 

53.9+3 

— 1 00 +O ' 13 
l.UU_o 13 

58.95 ±0.17 

1.13 


4 

0 44100 + °'° 0020 

U.44IUU_ 0 . 000 27 

351+I 1 

H 99C+0-061 

u.^^o_ Q 062 

6.68 ±0.19 

0.97 

2014 August 03 
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0 4412917 +0 ' 0000015 

i _ 0 0000033 

4-79+il 

— 20 [ 5+ 2 - 7 
ZU.O_ 20 .0 

311.93 ±0.08 

1.96 
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0 4412163 +0 ' 0000027 

U.^lZlO9 _ 0 QQQ0027 

8 - 1+1 

—4.641+3 

39.56 ±0.13 

1.32 

2014 September 01 

1 

0 4412914 +0 ' 0000008 
u.^izyr^_ 0 0000012 

4-45+j° 

-22.2+-J 

1188.89 ±0.08 

3.75 


2 

n 441 999Q+0.0000023 
U.441zzzo_ 0 0000019 

6-7+7 

—4.55+42 

139.57 ±0.14 

1.97 


3 

fl 441 S9n+0.000017 
U.441ozU_ 0 000017 

55.71+ 

-0-82lg.il 

39.22 ±0.17 

1.35 


4 

0 44050 +0 ‘° 0014 

U.44UOU_ q Q0017 

328+ 9 

H 9^9+0-049 

u.zoz_ 0 05 7 

17.96 ±0.19 

1.06 

Average 

1 

0 44129230 +0 ‘ 00000039 

U.^IZyZoU — 0.00000044 

4 96 +0 ‘ 15 
^• y °-0.16 

-20.38++ 

2417.76 ±0.08 
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0 4412209 + °' 0000011 
u.^izzuy_ 0 0000011 

‘ - DO -0.66 

_ 4 iq+0-27 
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3 

0 441819 +0 ' 000010 

U.'±'±iOi»_Q 000010 

54-21 l A 0 

_n 900 +O ' 074 

u.yuU-Q 075 

116.76 ±0.18 

1.68 


4 

0 44061 + °‘° 0014 

U.44UOI _ 0i00014 

3381 H 

—0 197+0- 030 
-0.031 

24.23 ±0.20 

1.30 
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Figure 5. The best-fitting models of the BETA spectra for each of the three observing epochs and the average spectrum. The radial velocity axis is 
given relative to the rest frame defined by the peak optical depth (vertical dotted line). The grey line represents the spectral data and the solid black 
line represents the best fitting line model. The dashed green lines denote the individual Gaussian components, which are identified by numbers that are 
ordered by descending peak line strength. The solid red line denotes the best fitting residual. 
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spectrum towards PKS B1740—517, averaged over all three 
epochs. It is clear that the quality of the spectrum in this band 
is extremely high and we did not need to excise any large re¬ 
gion due to RFI. As described above, only minor birdies were 
found, and these were due to hardware issues that either have 
been fixed or will not be present for the full ASKAP system. 
These observations make us optimistic that ASKAP will be a 
successful instrument in this frequency range. 


3 RESULTS 

3.1 HI absorption towards PKS B1740-517 

We search for the signature of HI absorption in the BETA data 
using a spectral line detection and modelling technique based on 
Bayesian model comparison (see |Allison et al.|2014| and refer¬ 
ences therein). This technique is a robust method of model se¬ 
lection, enabling us to both assign a significance to the detection 
and compare between increasingly complex model representa¬ 
tions of the data. The spectral data are modelled by convolving 
a physical model of the line, consisting of multiple Gaussian 
components, with the known gain response of the fine filterbank 
channels. Using this method we detect a single absorption line in 
the 21-cm spectrum of PKS B1740—517, located at an observed 
frequency of 985.5 MHz, equal to a HI redshift of z = 0.4413. 
By sequentially increasing the number of Gaussian components, 
and comparing the Bayes factors, we arrive at a best-fitting spec¬ 
tral model for each of our observing epochs, and for the average 
spectrum. In Table[2] we summarize the model parameters and 
in Fig.[5]we show the best fitting Gaussian components and the 
total model. 

We find that the combined data are best fitted by four com¬ 
ponents, corresponding to two narrow components (widths ~5 
and 8 km s -1 ) at the position of peak absorption and two broad 
components (widths ~50 and 350 km s -1 ), which are, respec¬ 
tively, redshifted and blueshifted with reference to the peak ab¬ 
sorption. The latter blueshifted broad component, with a peak 
at approximately 0.2 per cent of the continuum, is arguably the 
most tentative. However, we believe the veracity of this detec¬ 
tion for the following reasons: (a) the inclusion of this compo¬ 
nent is strongly warranted by the data above the formal noise, 
with an increase in log evidence of A\n(Z) = +24.23 ± 0.20 
and a decrease in reduced chi-squared for the best-fitting model 
parameters of Ax 2 /d.o.f. = —0.38; (b) other such features are 
not present elsewhere, either in this spectrum or other sources in 
the field (see Section [+3] >; (c) the feature is stable to changes in 
our continuum subtraction procedure. While formally detected 
above the noise, the peak strength of this feature is weak and 
further observations are needed for confirmation. We discuss the 
possible physical interpretations of these results in Section[5] 



Figure 6. The probability of detecting at least one intervening ab¬ 
sorber with column density greater than A+i in the PKS B1740—517 
field observed with BETA. The blue line represents a simple Gaus¬ 
sian model for the Hi gas, assuming that Cf = 1, a spin temperature 
of T S pi n = 100 K and Auso = 30kms -1 . The green line repre¬ 
sents the two-phase temperature-shielding model by |Braun| ( |2012) , with 
No = 1.25 x 10 20 cm -2 and Nqq = 7.5 x 10 21 cm -2 . The vertical 
grey line denotes a Damped Lya absorber (A+i > 2 x 10 20 cm -2 ), 
showing that with BETA we are only sensitive to systems with signifi¬ 
cantly higher column densities. 


radio source components with respect to the foreground ab¬ 
sorber ( |Briggs|1983|), int erstellar scintillation in the Milky Way 
Galaxy (Macquart 2005), and microlensing of the background 
source ( [Lewis & Ibata 2003). Due in part to limited channel 
sampling of the strongest absorption components, our data are 
sensitive (at 3 a) to fractional variations in the optical depth 
of only 30 per cent over a 2 month interval. For a background 
source traversing the absorber at a velocity close to the speed of 
light, our data are sensitive to fluctuations in the 21-cm opacity 
of more than 30 per cent on transverse scales less than 0.05 pc. 
If we consider instead the case of interstellar scintillation, for 
an absorber at z = 0.44, the required physical scales are much 
less than a few parsecs. Since we find no significant evidence 
for variability in the component parameters given in Table[2]this 
suggests that the bulk of the cold HI gas in these components 
is distributed uniformly over parsec scales, consistent with that 
found in Local Group galaxies (c.f. ~ 100 pc; |Braun|20 1 2) and 
radio galaxies ( Curran et al. 2013a). Given our time sampling of 
approximately a month, it is also possible that significant scin¬ 
tillation may occur on day and intraday time-scales to which our 
data are not yet sensitive. 


3.2 Spectral variability 

Our observations with BETA were carried out in three sepa¬ 
rate epochs with adjacent intervals of 40 and 29 d. Over these 
time intervals we can test for variability of the individual spec¬ 
tral components. Limited observational evidence in the literature 
indicates that ~ 10 per cent stochastic fluctuations may occur in 
the relative fluxes of some HI absorption components on time- 
scales of days to several weeks; prominent examples are the 
two intervening absorption systems towards the flat-spectrum 
quasar s AO 0235+164 ([Wolfe et al.|1982[ > and PKS B1127-145 
( [Kanekar & Chengalur||200l| , which exhibit flux variation but 
no significant shift in position or width. Several possible mod¬ 
els exist to explain such behaviour, which include motion of 


3.3 Intervening absorption towards other sources 

The field of view available with nine PAF beams enables a 
wider search for intervening Hi absorption in the spectra of 
multiple sources. Based on a nominal noise per channel of 
20 mJy beam -1 for the data averaged over all three observations, 
we searched the spectra of 72 sources brighter than lOOmJy 
within 1° of each beam centre, selected from the SUMSS 
( |Mauch et al.||2003| ) and MGPS-2 ( [Murphy et al.||2007] ) cata¬ 
logues, giving a 5 a detection limit for absorption against the 
weakest sources at the beam centre. Using the spectral-line find¬ 
ing technique of |Allison et al.|(2014| > we found no further Hi 
absorption in the spectra of radio sources in this field i.e. all 
spectra are consistent with the noise. 

We can estimate the expected number of intervening ab- 
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sorbers n with a column density greater than Nui by 

n(>N H i) = r f f(K 1 ,X)dXdK 1 
J n h1 J 

OO 

- EE / (^Vim 5 X) 5X 5Nfti , (3) 

Nm 

where / is the frequency distribution of systems with column 
density iVm and X is the comoving path length given (in a flat 
ACDM universe) by 

X{z ) = -Ev/fiMU + W 3 + «A. (4) 

To estimate / as a function of redshift we perform a simple lin¬ 
ear interpolation between distributions measured for the local 
Universe ( [Zwaan et al.|2005| ) and at 2 = 3 (Noterdaem e et al.| 
[2009] ). 

For each sight line towards our 72 sources we calculate the 
comoving path interval by 


jX(z ± 0.5 Sz) - X(z - 0.5 Sz), 

l°> 


if AW > Nsa, 
otherwise, 


(5) 


where Ns# is the column density sensitivity for detection of an 
absorption line as a function of spectral channel, and Sz is the 
channel separation in redshift. Given that the BETA spectra span 
a wide redshift interval, and that absorption can only be detected 
between the source and observer, we must also take into consid¬ 
eration the distribution of source redshifts. We therefore use the 
following redshift distribution model of |de Zotti et al.| ( [2010| ) 

A/'s(z) = 1.29 + 32.37z - 32.89z 2 + 11.13z 3 - 1.125z 4 , ( 6 ) 


which is determined from the population of radio sources 
brighter than lOmJy at 1.4 GHz. For any given sight line, the 
probability of the source being located at a redshift beyond z 
is equal to 


Pr(z 3 > z) 


J z °° A/' s (z')dz / 

Jo°° V s (z')dz'’ 


(7) 


and so the expected value for the absorption comoving path 
length probed by each spectral channel is given by 


(SX(z)) = Pr(z s > z)SX(z). 


( 8 ) 


We consider two models that convert the observed spectrum 
(in units of fractional absorption) to the column density sensitiv¬ 
ity. First we use a simple conversion assuming line-of-sight gas 
of spin temperature T sp i n = 100 K, covering factor Cf = 1 and 
velocity FWHM Auso = 30kms _1 . The column density sensi¬ 
tivity is then given by 


N 5<7 = 5.82 x 10 21 


' Tspin ‘ 


Au50 

100 K 


30 kms -1 


(9) 


where is the optical depth sensitivity estimated from the 
spectral noise per channel and assuming a covering factor of 
unity. In the case of our second model we use a simple two-phase 
sandwich geometry for the gas where the cold neutral medium 
(CNM; T c ~ 100 K) is sandwiched by a layer of warm neutral 
medium (WNM; T w > 5000 K), which acts as a shield from lo¬ 
cal high energy UV and X-ray radiation (for further details see 
|Kanekar et al.|20ll| and |Braun 12012 ). In this case the column 
density sensitivity is given by 


N 5a = No + 




( 10 ) 


where N 0 is the threshold column density of WNM at which 
the CNM can form, and Noo is the saturation column density 


Table 3. Summary of emission line properties from our Gemini-South 
observations. 


Line log(flux) Equivalent width log(luminosity) 

(ergs -1 cm -2 ) (A) (ergs -1 ) 


[On] A3727 

-16 11+ 004 

23.60 ±2.26 

40.74t°° 4 

H (3 


4.22 ± 1.00 

40.521°'°® 

[O ill] A4959 

— 16 nn+ 0 - 04 
io.uu_ 0 05 

8.52 ±0.91 

40.86l°° 4 

[O hi] A5007 

—15-58t°-°l 

23.73 ±0.95 

41 27 +0 ' 02 
^• Z '-0.02 

[Oi] A6300 

—16.07±°;“ 

6.68 ±0.96 

40.78lo.o7 


as the optical depth in CNM tends to infinity. We use values 
for these parameters of N 0 = 1 .25 x 10 20 cm -2 and Noo = 
7.5 x 10 21 cm -2 , obtained by Braun ([2012]) from model fit¬ 


ting to detailed 21-cm observations of the Milky Way Galaxy, 
Messier 31 and the Large Magellanic Cloud. 

The probability of detecting J\f absorbers with column den¬ 
sities greater than Nm is then given by 




Pr[A/"(> Nhi)] = f- exp(- 


-n), 


(ID 


where n is the expected number of absorbers given by Equa¬ 
tion]^ Based on the above models of the Hi gas we expect to 
detect less than ~ 0.05 intervening absorption systems with col¬ 
umn densities greater than about 10 21 cm -1 , with a correspond¬ 
ing probability of less than ~ 0.05 for detecting at least one such 
system (see Fig.[ 6 ]). Therefore the lack of any intervening detec¬ 
tion in our data is consistent with the known distributions of HI 
in the local and distant Universe. 


4 MULTI-WAVELENGTH FOLLOW-UP 

4.1 Optical spectroscopy with Gemini-South 

4.1.1 Observations and data analysis 

As discussed in Section [±2] prior to our observations with 
the BETA telescope the redshift of the optical counterpart to 
PKS B1740—517 was not well known, with no spectroscopically 
determined value available from the literature. Furthermore, the 
HI redshift did not match any of the published photometric esti¬ 
mates. The interpretation of the observed absorption required an 
accurate knowledge of the host galaxy redshift, so that the HI 
gas could be associated either with the host galaxy of the radio 
source or with an intervening system. 

We therefore obtained time from Gemini-South under the 
Director’s Discretionary Time program (proposal code GS- 
2014B-DD-2) to acquire long-slit spectroscopy and gW-band 
imaging with the Gemini Multi-Object Spectrograph (GMOS; 
[Hook et al.|2003| ). The long-slit spectroscopy was taken with the 
R400 grating, utilizing two grating settings with central wave¬ 
lengths 7000 and 7050 A. This was chosen to put the U/3 and 
[O III] lines in the centre of the spectrum for a host redshift equal 
to that of the HI absorption. At each grating setting we obtained 
3 x 900 s exposures, giving a total exposure time of 1.5 h. A 
1.5-arcsec slit was used, and the CCD was binned 2x2 on-chip. 

The spectra were reduced in I RAF, using standard tech¬ 
niques within the Gemini package. The spectra were bias- 
subtracted and flat-fielded using calibration frames from the 
Gemini Facility Calibration Unit (GCAL). Wavelength calibra¬ 
tion made use of Cu-Ar comparison spectra, and the spectrum 
was flux calibrated using observations of the standard star Feige 
110 . We estimate an rms error of 0.9 A in the wavelength calibra¬ 
tion, based on the spread in the Cu-Ar lines. We also acquired 
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Figure 7. A three-colour optical image constructed from gW-band observations with the 8-m Gemini-South telescope. PKS B1740—517 is identified 
with the red galaxy located between the two solid horizontal bars. The orientation and position of the 1.5-arcsec slit are indicated by the parallel lines. 
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Figure 8. The optical spectrum towards PKS B1740—517 from long-slit observations using the GMOS instrument on the Gemini-South Telescope. The 
black line denotes the full spectrum and the vertical dotted lines denote the rest wavelengths of known species. We have performed a model fit to the 
stellar population (red line) using pPXF, where the emission lines and Telluric regions (blue lines) have been masked, obtaining a systemic redshift of 
Zsys = 0.44230 ± 0.00022. 
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log([0ll]X3727/H(S) 

Figure 9. The blue classification diagram for redshifted emission line 
galaxies |Lamareille|2010| , where the position of PKSB1740—517 is 
indicated by the error bars. 


images with GMOS in the g'r'i bands. Each image consisted 
of 4 x 75 s exposures, which were binned 2x2 on-chip, giving 
a pixel scale of 0.16 arc sec pixel -1 . A three-colour image ( g' 
as blue, r' as green and i! as red) is shown in Fig.[7] Adjusting 
the original astrometry of the Gemini image using the B } -band 
image from the SuperCosmos Sky Survey ( [Hambly et al.|200l| 
we find that the centroid position of the large central galaxy (in¬ 
dicated by the solid horizontal bars in Fig. [7} agrees with the 
ICRF2 position of PKS B1740—517 to within 0.6 arcsec. 

The full GMOS spectrum is shown in Fig.[8]with a mea¬ 
sured FWHM spectral resolution of 11 A from fitting the Cu-Ar 
arc spectra. Clearly apparent are the [O III] AA4959, 5007 dou¬ 
blet, the [O II] A3727 and [O I] A6300 emission lines, all consis¬ 
tent with a redshift that matches the HI absorption. There is also 
evidence of [O III] A4363 emission. However, it is contaminated 
by subtraction of the 6300 A sky line and so we cannot infer 
accurate kinematics for this line. We note that there is almost 
no sign of Baimer-line emission, with weak R/3 at 4861A and 
Ha redshifted out of the observed band. In order to obtain an 
estimate of the systemic redshift of the host galaxy, we use the 
Penalized Pixel-Fitting method (pPXF;|Cappellari & Emsellem 


obtain a systemic redshift of 2 = 0.44230 d= 0.00022, where the 
uncertainty is derived by summing in quadrature the formal er¬ 
ror due to the noise (cr Zj noise = 0.00013) and the wavelength 
calibration (a ZjCa 1 = 0.00018). We therefore conclude that the 
HI absorption system seen towards PKS B1740—517 is intrinsic 
to the galaxy hosting the radio AGN. 


2004 to fit stellar population synthesis models ( |Vazdekis et al. 
2010) to the underlying stellar continuum. Using this method we 


4.1.2 Emission line ratios and kinematics 

In Table[3] we summarize the properties of individual emis¬ 
sion lines seen in the Gemini spectrum, where the fluxes and 
equivalent widths are estimated from fitting Gaussian models 
to the emission lines after subtraction of the best-fitting stel¬ 
lar population models. The continuum level, used in calculating 
the equivalent width, is measured adjacent to the relevant line 
from the original unsubtracted spectrum. Equivalent width ra¬ 
tios of [O iii]A 5007/H/3 = 5.60 ± 1.43 and [O n]A3727/H/3 = 
5.64 d= 1.34 place PKS B1740—517 in the Seyfert region of the 
blue classification diagram for redshifted emission line galax¬ 
ies (Fig.[9j |Lamareille||2010| >. The [OI] A6300 line is particu¬ 
larly strong ([O l]A6300/[0 III] A5007 = 0.28 d= 0.04), which in 
combination with the parity in strength of the [O III] A5007 and 



-2000 -1500 -1000 -500 0 500 1000 1500 2000 

Velocity offset (km/s) 


Figure 10. The HI 21-cm absorption, oxygen and H/3 emission lines in 
the rest frame defined by the systemic redshift of the galaxy ( z sys = 
0.44230 =b 0.00022). In all cases the continuum has been subtracted 
from the spectrum. The vertical dashed lines denote the location of 
HI absorption components 1 (left) and 3 (right), which correspond to 
the most prominent features. The red horizontal error bar indicates the 
=L46kms _1 uncertainty in the systemic velocity. 


[OII] A3727 lines could be an indicator of Low-Ionization Nu¬ 
clear Emission line Region (LINER; |Heckman| 1980| ) behaviour. 
Strong [O I] A6300 and [O III] A4363 emission arise from exci¬ 
tation by fast shocks (300 - 500 km s -1 ), which also indicates 
possible interaction between the radio jets and interstellar gas 
( [Dopita & Sutherland! 1 995). 

In Fig.|10|we show~the velocity structure of the optical 
emission lines (after subtraction of the fitted stellar component) 
in the rest frame defined by the systemic redshift. Despite the 
low spectral resolution (ranging from 350 to 650 km s -1 across 
the observed band), there is strong evidence for double-peaked 
structure in the [O III] and [O I] A6300 emission lines. In particu¬ 
lar, the high S/N [O III] A5007 line clearly shows two peaks sep¬ 
arated by a velocity of approximately 500 km s -1 . Conversely, 
the strong [OII] A3727 line appears to have only a single peak 
at the systemic velocity, but is at much lower S/N than the two 
[O III] lines and is likely confused by the unresolved doublet. 

The relative strength of the narrow oxygen and R/3 lines 
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Spectral pixels (0.148 nm/pix) 
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Figure 11. Top: an extract of the two-dimensional spectrum showing just 
the [O in] A5007 line. The axes show the pixels from the arbitrary start 
of the cut-out. The spatial axis pixels increase in a roughly northerly 
direction (consult Fig .[7] for the orientation of the slit on the sky). The 
two apertures used to examine the spatial profile of the emission line are 
shown, one for each peak of the A5007 line. Bottom: the spatial profile 
averaged over each aperture, after subtracting the average profile across 
the entire cut-out. The colour of the profile matches the aperture shown 
in the 2D spectrum. Also plotted (in black) are Gaussian fits to the profile 
used to estimate the peak location for each pixel. The profiles are plotted 
as a function of distance from the peak of the average profile, with the 
pixel values indicated for comparison (positive offsets are in a roughly 
northerly direction). 


strongly indicate that the emission is dominated by gas ionized 
by the AGN, rather than star formation. So far only a few AGN 
with double-peaked [O III] lines have been reported in the liter¬ 
ature ( |Liu et al.|2010| [Smith et al.|2010||Shen et al.|2011| ), and 
their origin has largely been ascribed to kinematic effects within 
a single, rather than binary, AGN (see |Fu et al.|2012| . In Fig.pT] 
we show the spatial distribution of the [O III] A5007 line along 
the slit, which exhibits a clear separation of 0.22 arcsec (1.3 kpc 
at z = 0.44) between the peaks. Such behaviour could be due 
to an outflow of AGN-driven gas (e.g |Fischer et al.|2011| ), per¬ 
haps even directly as a result of radio jet-cloud interaction (e.g. 
|Stockton et al. 12007] |Rosario et al.|2010) , but could also be due 
to kpc-sc ale rotation of the ionized gas. While detection of an 
ionized outflow would be typical of the host galaxy of a power¬ 
ful compact radio source, the large scale seen here would be at 
odds with other GPS sources, where ionized outflows are seen 
on the same scale as the radio jet (e.g. |Holt et al.|2009) . Several 
recent integral field spectroscopic (IFS) observations of lumi¬ 
nous (radio-quiet) Type II AGN (e.g. |Liu et al.|2013| |Harrison| 
|et al.|2014]|McElroy et al.|2015| ) have confirmed the prevalence 
of large-scale AGN-driven winds. A similar spatially resolved 


IFS study of the narrow line emission would confirm the origin 
of the ionized gas kinematics in PKS B1740—517. 


4.2 Archival data from XMM-Newton 

4.2.1 Data analysis 

By matching the ICRF2 position of PKS B1740—517 to 
sources listed in the 3XMM-DR4 XMM-Newton Serendip¬ 
itous Source Catalogue, we found a bright X-ray source, 
3XMM J174425.3—514444, at an angular separation of 1 arcsec. 
This field was serendipitously observed as part of targeted obser¬ 
vations towards HD 160691, a fifth-magnitude star also known 
as /xArae, which is 5.9arcmin away from PKS B1740—517. 
Via the XMM-Newton Science Archiv^] we obtained four X- 
ray spectra of PKS B1740—517 (see Fig.[l2]), from two sepa¬ 
rate observations of integration time 7800 and 11 700 s, and fit¬ 
ted these simultaneously using the modelling package XSPEC 
( |Arnaud||1996| ). We found that the X-ray spectra are well fit¬ 
ted by a standard absorbed power-law model. Because of the 
relatively low Galactic latitude of b = —11? 5, we incorpo¬ 
rate an additional fixed absorption component in our model 
due to the Galactic foreground, estimated within XSPEC to 
be A^H,Gai = 1-09 x 10 21 cm -2 . Using this model, we 
obtain an estimate of the intrinsic column density towards 
PKSB1740-517 of N u ,x = 1.211^3 x 10 22 cm -2 , a pho¬ 
ton index of T = 0.801°'24 and a 1 keV normalization of 
4-8ti.5 x 10 -5 photons cm -2 s -1 keV -1 . While the intrinsic 
column density is consistent with that measured for other X- 
ray observed GPS radio sources, the photon index is hardej^] 
possibly indicating an obscured Compton-thick AGN, similar 
to that seen for PKS B1607+26 fTengstrand et al.|2009| . If the 
observed spectrum is indeed arising from reflection of the pri¬ 
mary X-ray emission off a Compton-thick obscuring medium, 
then sight-lines to the AGN may instead have column densi¬ 
ties in excess of 10 24 cm -2 . Using the above best fitting model 
to the XMM-Newton spectra we estimate an X-ray luminosity 
of L 2 - iokeV ~ 4.3 x 10 43 ergs _1 , consistent with the dis¬ 
tribution of luminosities measured for other GPS radio sources 
( |Tengstrand et al.|2009| . 


4.2.2 An X-ray upper limit on the HI spin temperature 

Combining the X-ray derived total column density of elemental 
hydrogen and the 21-cm opacity data allows us to estimate an 
upper limit on the HI spin temperature, and compare this with 
direct measurements of the spin temperature in other galaxies. 
In the absence of any Hi emission, the 21-cm optical depth is 
given by 

r 2 i = -ln(l + -|^A, (12) 

where A S' is the change in flux density due to absorption, S CO nt 
is the continuum flux density and a is the covering factor, i.e. 
the fraction of background source obscured by the foreground 
absorber. The column density of HI gas (in units of atoms cm -2 ) 
is then given by 

N m = 1.823 x 10 18 T spin J r 21 (v)dv, (13) 

where T sp i n is the 21-cm spin temperature (in units of K) and 
/ T21 (v)dv is the rest-frame velocity integrated optical depth (in 


7 http://xmm.esac.esa.int/xsa/ 

8 (F) = 1.66, or = 0.36 for the sample of 16 X-ray observed GPS 
sources discussed by jTengstrand et al.]{2009| . 
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Figure 12. Upper panel: four X-ray spectra towards PKS B1740—517, 
observed with the XMM-Newton satellite. The solid lines represent the 
best-fitting model produced using an absorbed power-law in XSPEC. 
Lower panel: the residuals from the best-fitting models, indicating scat¬ 
ter consistent with the noise. 


units of kms -1 ). For sight lines where several gas components 
may exist at different spin temperatures, T sp i n is the column- 
density-weighted harmonic mean over those components. As¬ 
suming that the 21-cm and X-ray spectra probe similar sight 
lines through the gas, an upper limit on the spin temperature can 
be obtained from 


T S pi n < 5485 


iV H ,x 

f T2l(v)dv~ 

10 22 cm -2 

1 kms -1 


K, 


(14) 


where Nn,x is the total column density of elemental hydro¬ 
gen (in atomic, ionized and molecular form) estimated from the 
XMM-Newton spectra. 

We estimate the 21-cm optical depth from our BETA data 
using Equation|T2] which of course depends on the value we 
adopt for the covering factor (cf). Assuming a covering fac¬ 
tor of unity (which gives a lower limit to the optical depth and 
hence a spin temperature upper limit), our best estimate of the 
total integrated optical depth from the average BETA spectrum 
is2.73±0.10 km s -1 . Based on the column density of total hy¬ 
drogen derived from a simple absorbed power-law model of the 
X-ray spectra, we obtain an upper limit on the spin temperature 
of T S pin < 2430 =L 1050 K. The large fractional uncertainty in 
the X-ray derived Ah,x results in a poor constraint on this up¬ 
per limit, yet is consistent with other direct measurements of ex- 
tragalactic spin temperatures, both from 21-cm observations of 
known damped Lya absorption systems towards quasars (Curran| 
|2012| [Kanekar e tal.|2014| ) and by combining the intrinsic HI ab¬ 
sorption in quasars with extinction estimates from their optical 
continuum (e.g. PKS B1549—790; [Holt et al.[2006| . While this 
agreement is reassuring, the assumption of similar sight lines to 
the radio and X-ray sources is tentative at best, depending on the 
extent of the radio source and mechanism generating the X-ray 
emission. 


5 DISCUSSION 

The HI absorption line profiles associated with radio galaxies 
are typically complex, encoding small-scale variation in both the 
21-cm opacity and the background source structure that are often 
difficult to disentangle and accurately model. However, we can 
attempt to overcome this by amalgamating information from the 


integrated HI spectrum, the radio continuum and any available 
multiwavelength data. In the following we discuss the inferred 
properties of the background radio source, followed by the in¬ 
dividual components seen in the absorption spectrum and their 
most likely physical interpretation. 


5.1 Properties of the background radio source 

PKS B1740—517 is a powerful young GPS radio source, which 
through optical imaging and spectroscopy we have established 
is hosted by a luminous Seyfert galaxy at z = 0.44. A least- 
squares fit to the spectral energy distribution (SED) shown in 
Fig.|T] yields a spectral age of approximately t age « 2500 yr, 
where we have assumed an equipartition magnetic field ( |Murgia| 
|2003| . From VLBI observations, King found that the 2.3-GHz 
flux density is contained within two 10 mas-scale radio com¬ 
ponents, in a ratio of approximately of 3.9 to 1 fKmg][T994) . 
The angular separation of the two components is 52 mas, corre¬ 
sponding to a physical scale of ~ 300 pc. From the integrated 
SED alone we cannot confirm whether the components repre¬ 
sent hot-spots associated with two jets or a core and jet. In pre¬ 
vious VLBI studies of compact sources ( [Stanghellini et al.|1997[ 
|Snellen et al.|2000||de Vries et al.|2009) it was found that GPS 

radio sources associated with galaxies (as opposed to quasars) 
tend to exhibit compact double (CD) or compact symmetric ob¬ 
ject (CSO) morphologies. In the following discussion we there¬ 
fore assume that the VLBI-scale structure seen in this source 
represents radio emission associated with jet hotspots, but can¬ 
not rule out the possibility of a core and jet without further spec¬ 
tral information for the two source components. 

The brightness asymmetry of the two jet components could 
be either caused by an intrinsic difference in the radio luminos¬ 
ity, possibly resulting from an inhomogeneous circumnuclear 
medium, or by an orientation effect. In the latter case we can 
estimate the angle i at which the radio jet axis is inclined to the 
line-of-sight, as follows 


7 2 

1/2 

- R l/(2-a) 

~ f 

7 2 — 1 


R l/(2-a) 

FI, 


where Rj c is the brightness ratio, a is the spectral index and 
7 is the bulk Lorentz factor in the jet. Assuming that the in¬ 
trinsic luminosities of the two jet components are equal and 
that a = — 0.74, 77j c = 3.9 (and a typical Lorentz factor in 
the range 5 < 7 < 10), we estimate that the jet axis is in¬ 
clined by i ~ 76° with respect to the line-of-sight, giving a 
jet axis radius of r ~ 150 pc. Based on the spectral age we 
estimate an average separation velocity of approximately 0.4 c 
over the lifetime of the radio AGN, consistent with directly mea¬ 
sured expansion velocities of similarly powerful GPS and CSS 
sources (see [de Vries et al.||2009| and references therein). The 
combined linear size, X-ray luminosity and 5-GHz radio lumi- 
nosit 300* GHz ~ 8.7 x 10 43 ergs 1 ) are all typical of other 
GPS radio sources in the literature (e.g. |Tengstrand et al.|2009| . 


5.2 The cold neutral gas in PKS B1740-517 

5.2.1 Absorption components 1 and 2 

Absorption components 1 and 2 in the BETA spectrum have sim¬ 
ilar velocities and may either be separate cold clouds of HI or 
the same structure seen against both source components. With 


9 For comparison with the literature we define the 5-GHz luminosity as 
Lb GHz = GHz £5 GHz, where i/ 5 GHz is the rest-frame frequency 
and S 5 ghz is the luminosity density. 
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a width of Avr.u < 5kms _1 , the deep component is unusually 
narrow for extragalactic HI absorption associated with the host 
galaxies of radio AGISrj Typical widths often range from more 
than 10 to 1000 km s^^and are largely attributed to turbulent 
and bulk motion of gas that is either rotating or radially outflow¬ 
ing/infalling with respect to the nucleus (see e.g. | Allison et al.| 
2013 Gereb et al. 2015 and references therein). Narrow absorp¬ 
tion is often found near to the systemic velocity and thought to 
be associated with gas in a regularly rotating structure, such as a 
disc or ring, through which a single sight line to a compact radio 
source picks out the velocity dispersion. 

In Fig.[K)]we compare the velocity structure of the Hi ab¬ 
sorption and optical emission lines. The prominent narrow com¬ 
ponent is blue shifted by approximately —200 d= 50 km s -1 with 
respect to the systemic velocity and also the dynamical centre 
of the [O III] emission lines. Assuming that the radio source 
is co-spatial with the AGN, this offset would indicate that the 
gas is either rotating or outflowing radially towards the observer 
from the nucleus. A similar 200km s -1 velocity offset is seen 
for narrow Hi absorption (Auso ~ 90kms _1 ) in the nearby 
Seyfert 2 galaxy NGC2110. Detailed modelling by |Gallimore| 
|etal.| ( [1999| ) showed that the absorbing gas is located in a ~ 200- 
pc disc, close to the radio source and with similar orientation to 
the jet axis. The velocity offset seen here would also suggest that 
the separation of the gas and radio source is small, indicating a 
location that is ~ 150 pc away from the central nucleus. Since 
a jet-driven radial outflow would likely result in significant tur¬ 
bulence of the gas, and therefore the narrow absorption would 
represent an intact cloud that has become entrained, it is more 
likely that we are seeing absorption through a foreground struc¬ 
ture. 

It is also possible that an unforeseen systematic error in 
the Gemini spectrum (e.g. contamination from dust obscuration) 
may have caused an artificial offset between the optical redshift 
and the narrow absorption. If the narrow feature is instead at the 
systemic velocity, this would place the absorbing gas at a dis¬ 
tance much further out in the galaxy (similar to PKS 1814—637; 
|Morganti et al.|20lT) , and would represent tangentially rotating 
gas seen in absorption against the central radio AGN. Similarity 
between the line ratios of absorption components 1 and 2 and 
the flux densities of the radio components would also suggest 
lines of sight to each. Under this scenario the distance of the 
absorbing gas from the nucleus would simply be given by 


^disc 

Au a bs 


1 4kpc, 


(16) 


where r source is the physical separation of the two radio com¬ 
ponents, Udisc is the rotational velocity of the gas (assumed to 
be ~ 200kms _1 ) and Au a bs is the velocity separation of the 
two absorption components. However, there is no obvious ev¬ 
idence of the effect of dust obscuration in the Gemini data and 
verifying this would require deeper optical imaging and resolved 
spectroscopy. 

The peak HI component is sufficiently narrow that we can 
estimate a reasonable upper limit on the kinetic temperature of 
the gas. We use the following relationship between the tempera¬ 
ture and line width 

Tl ~8E(2)l| 1(a,, »" )> - 538±35K ’ 071 

where Auso is the FWHM of the spectral line, mu is the mass of 
a hydrogen atom and ku is the Boltzmann constant. Given that 
turbulent motion of the gas will result in line broadening, we 


expect that the true kinetic temperature could be significantly 
less than this and in a phase similar to the CNM of the Milky 
Way (Tk ~ 40-200 K; |Wolfire et al.|2003) . Assuming that the 
gas is dense enough to be collisionally excited, the kinetic tem¬ 
perature is a reasonable proxy for the HI spin temperature, and 
provides a stronger constraint than the X-ray upper limit derived 
in Section [4,2.2| However, irradiation by UV emission from the 
AGN and continuum emission from the powerful radio source 
may drive the spin temperature to much higher values ( jBahcall| 
|& Ekers|1969| >. Given the high luminosity of the radio source, 
and the fact that we detect HI absorption, we find it improba¬ 
ble that the spin temperature could be dominated by 21-cm ra¬ 
diation. We therefore assume that the gas is sufficiently dense 
and/or shielded from the continuum to be collisionally excited, 
but acknowledge that otherwise the spin temperature could be 
an order of magnitude higher than the kinetic temperature. 

If the peak absorption is located in front of the fainter radio 
component it would imply that all of the radio flux is being ab¬ 
sorbed (given the flux ratio of 3.9 to 1 at 2.3 GHz), in which 
case we would expect the profile to be saturated. Given that 
the ASKAP spectral channels do not fully sample this narrow 
lin£3 we cannot adequately resolve the profile to test for sat¬ 
uration. However, resolved HI absorption lines associated with 
radio AGN do not typically exhibit saturation and such a sce¬ 
nario would indicate that the gas here is considerably colder than 
the upper limit given in EquationffT] We therefore favour the 
scenario that the peak absorption is seen in front of the brighter 
radio source. In which case the HI covering factor of the total ra¬ 
dio flux is Cf = 0.8, and so the velocity integrated optical depth 
for the deep component is 1.55 ± 0.07 km s -1 . Using the above 
upper limit to the kinetic temperature, and the relationship given 
by EquationpG] we estimate a corresponding column density of 
Nui < (1.52 =b 0.12) x 10 21 cm -2 . The angular size of the 
brighter radio component measured by |King| < |1994| ) is 6.7 mas 
x 5.8 mas (at 2.3 GHz), which at a redshift of z = 0.44 equates 
to a physical scale of 38 x 33 pc 2 . We therefore estimate a mass 
of Mm ~ 5 x 10 4 M© for the foreground gas cloud. Uncertain¬ 
ties in the spin temperature and covering factor mean that the 
column density and mass could be an order of magnitude larger. 


5.2.2 Absorption components 3 and 4 

Component 3 (Fig.[5j is consistent with the dynamical centre 
of the [O ill] emission and offset by approximately —100 zb 
50kms _1 with respect to the systemic velocity. Its width 
(AU50 ~ 50kms _1 ) is comparable with Hi absorption often 
seen through the inclined discs of nearby Seyfert galaxies (e.g. 
|Dickey||1982| |Gallimore et al.||1999| | Allison et al.||2014| >. It is 

likely that component 3 is arising from tangentially rotating HI 
gas in a disc or ring, which would give a radial velocity close to 
that of the system. On the other side of the absorption complex 
we tentatively see a weak, broad 300km s -1 feature (compo¬ 
nent 4) shifted by —350 zb 50 km s -1 with respect to the system. 
Such a profile is typical of jet-driven neutral outflows seen in the 
host galaxies of other powerful radio AGN and we discuss this 
further in Section l53l 

Under the alternative scenario discussed above, whereby 
the deep narrow absorption is actually at the systemic redshift 
(and an unknown systematic error is causing a velocity offset), 
we infer a similar model to that of PKS B1814—637 ( [Morganti] 
|et al.pOTT) , where the narrow absorption is arising from gas 
in the larger galaxy and the broader features are the result of a 


10 |Borthakur et ab] j2014) also found very narrow intervening absorp¬ 
tion through the dwarf galaxy UCG7408 (Auso = 1.1 km s -1 ), corre¬ 
sponding to cold gas at a kinetic temperature of Tk ~ 26 K. 


11 Making this system an excellent future target for proposed high spec¬ 
tral resolution zoom modes with the SKA pathfinders. 


© 2015 The Authors, MNRAS 000,[l]{T8] 
































HI gas in a young radio galaxy using ASKAP 15 


rotating circumnuclear disc close to the radio source. However, 
there is no evidence so far to suggest that the optical redshift is 
being corrupted by a systematic error, and further observations 
at higher spatial and spectral resolution would shed light on this 
possibility. 


5.3 A possible jet-driven neutral gas outflow 


Jet-driven neutral gas outflows are direct indicators of the feed¬ 
back between a powerful radio AGN and the surrounding ISM 
( |Hardcastle et al.||2007| ). Gas outflows are found to be ubiqui¬ 
tous in the host galaxies of powerful compact radio AGN, and 
signify the growth of the source as it expands through a natal 
cocoon of dust and gas (e.g. |Holt et al.||2009) . In the case of 
PKSB1740—517, absorption component 4 (Fig.[5j may repre¬ 
sent a neutral gas outflow that is being driven by the approach¬ 
ing jet as it drives into the circumnuclear gas. The mass outflow 
rate of the neutral gas can be estimated by assuming a simple 
superwind model ( |Heckman et al.|2000} : 


Mhi — 30 


Q r Nm v 

47 t 1 kpc 10 21 cm -2 300 km s -1 


M© yr \ 


(18) 


where Q is the solid angle into which the gas is flowing (as¬ 
sumed to be ~ 7r), r is the radius, v is the velocity of the out¬ 
flow, and Nm is the column density of HI gas. Using an outflow 
velocity of v ~ 300kms _1 and a column density estimate of 
Nm ~ 1 x 10 21 cm -2 , which assumes that the outflowing gas 
has a similar spin temperature to that of the narrow absorption, 
we estimate an outflow rate of ~lM©yr _1 . This is strongly 
dependent upon our assumptions about the geometry and spin 
temperature of the gas, and we note that the outflow rate may be 
an order of magnitude higher. Within this framework the narrow 
absorption may either represent unperturbed circumnuclear gas 
or a dense cloud that has become entrained in the outflow. 


6 CONCLUSION AND FUTURE PROSPECTS 

Using commissioning data from the six-antenna BETA of 
ASKAP we have discovered Hi 21-cm absorption that is in¬ 
trinsic to the host galaxy of the powerful compact radio source 
PKSB 1740—517. Our result demonstrates the excellent poten¬ 
tial of ASKAP to search for HI in a continuous redshift range 
of z = 0.4 - 1.0, equivalent to look back times in the range 4.2 
- 7.7 Gyr. Such a capability is made possible by the excellent 
radio quiet environment of the observatory site at these frequen¬ 
cies. Furthermore, the detection and characterization of a com¬ 
plex absorption profile is very encouraging for the study of the 
neutral environments of radio AGN in future wide-field, broad¬ 
band absorption surveys with ASKAP and other SKA pathfind¬ 
ers. 

Optical spectroscopy with the 8-m Gemini-South telescope 
confirms that the HI absorption is intrinsic to the host galaxy of 
the radio source. The galaxy has a luminous active nucleus, with 
several forbidden oxygen lines exhibiting double-peaked struc¬ 
ture on large-scales that indicate possible AGN-driven outflows 
of ionized gas. Archival data from the XMM-Newton satellite 
show that the soft X-ray emission is absorbed by a dense obscur¬ 
ing medium, possibly Compton thick. The profile of the HI line 
is complex and exhibits the typical narrow and broad features 
associated with absorption in the host galaxy of a compact radio 
AGN. Based on the significantly blueshifted broad absorption 
we conclude that PKSB 1740—517 represents a recently trig¬ 
gered radio source ( t age « 2500 yr) in which the jets are break¬ 
ing through the surrounding cocoon of dense circumnuclear gas. 


This object requires further study, with spatially resolved opti¬ 
cal and 21-cm spectroscopy, and deeper optical imaging, to fully 
understand the neutral and ionized gas kinematics. 
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